Space-charge modulation in vacuum microdiodes at THz frequencies 
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We investigate the dynamics of a space-charge limited, photoinjected, electron beam in a micro- 
scopic vacuum diode. Due to the small nature of the system it is possible to conduct high-resolution 
simulations where the number of simulated particles is equal to the number of electrons within the 
system. In a series of simulations of molecular dynamics type, where electrons are treated as point- 
charges, we address and analyze space-charge effects in a /xm-scale vacuum diode. We have been able 
to reproduce breakup of a single pulse injected with a current density beyond the Child-Langmuir 
limit, and we find that continuous injection of current into the diode gap results in a well defined 
train of electron bunches corresponding to THz frequency. A simple analytical explanation of this 
behavior is given. 
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Vacuum microelectronic devices have shown some 
promise as sources of radiation with frequency ranging 
from GHz to THz. Examples of this include microtriodes 
PJ [2] , nanoklystrons [3] and traveling wave tubes [4] . 
Most such devices rely upon field emission as a source of 
electrons, although photoemission could also have similar 
advantages in terms of cold-cathode behavior. In cither 
case, it is reasonable to revisit some of the basic physics 
that describe macroscopic vacuum electronic devices and 
to see how they might manifest differently on the mi- 
croscale. One example is the problem of space-charge 
limited flow in a diode. The classic Child-Langmuir equa- 
tion 51 E] describes this to an extent. It is derived as- 
suming continuous flow of electrons from the cathode to 
the anode in a planar diode of infinite width. The initial 
analysis has since been improved upon to include the ef- 
fects of limited emission area [7J [S] finite pulse length 
[9], and quantum effects in small diode gaps [TO], to 
name a few. These improvements have been concerned 
with finding the onset of virtual cathode formation, so 
as to determine the upper limit of current flow that can 
be smoothly transmitted through a vacuum diode. In 
this study we look at the dynamics of the beam once the 
space-charge limit has been attained. 

The model used consists of a planar cathode with a cir- 
cular emitter area of radius of 250 nm seperated from a 
planar anode of infinite width by a vacuum gap of length 
d. Due to the small size of the system under study, it is 
possible to treat all electrons within the system individ- 
ually. An external static electric potential V is applied 
across the gap an gives rise to an electric field, E z = V/d, 
which accelerates the electrons from the cathode towards 
the anode. In addition to the external field an electron, 
i, also experiences the internal Coulomb field due to the 
space-charge created by the other electrons: 
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FIG. 1: Schematic model of the system simulated. The vari- 
able parameters are the applied voltage, V , and the gap width, 
d. The resulting currents at the cathode (emitter), jc{i), and 
at the anode (absorber), jA{t), are measured. 



where eo is the permittivity of free space, e is the electron 
charge, and the sum runs over all other electrons present 
in the gap. is the relative distance between electron i 
and j. Yielding the total field E tot = E sc + [0, 0, E x ] 

Dynamic behavior is obtained by solving the classical 
(Newton) equations of motion using a numeric Velocity 
Verlet integrator with a discrete time step of 1 fs. During 
a simulation electrons continuously enter and leave the 
simulation box, such that the number electrons in the 
device is not conserved. Relativistic effects and radiation 
caused by acceleration of electrons are safely neglected as 
all occuring velocities are much smaller than the speed 
of light. The model is shown schematically in Fig. [T] 
Note that the z-coordinate is defined by the normal to 
the surface of the cathode. 

The simulations conducted are divided into two groups 
according to the time profile of the laser pulse exciting 
photoclectrons. (I) Pulse scenario, in which an electron 
pulse of duration less than the transit time is created by 
varying the laser intensity for a fixed applied potential, 
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FIG. 2: Shape of the electron pulse when the externally 
applied voltage is varied. A fixed Gaussian photon profile 
is used. Only the highest voltage gives a 100% yield. In- 
set is earlier published experimental results, reproduced from 
Ref. [9j- The experimentally obtained double peaked shape 
was assigned to space-charge effects, this structure is recov- 
ered for the Medium voltage. 



or varying the the external potential for a fixed laser in- 
tensity. (II) Beam scenario, in which a continuous laser 
beam is emitted and either the gap length, or the external 
potential is varied. Reported values are from represen- 
tative simulations for the specific parameter-set. How- 
ever, all simulations have been repeated 10 times and 
only small variations are observed. 

The emission model used in the simulations assumes 
that photoelectrons are available in proportion to the 
photon intensity. Whether these photoelectrons are emit- 
ted is entirely dependent on whether the local electric 
field at their place of origin accelerates them away from 
the cathode or not. For an electron to enter the system 
at the cathode it must experience a force towards the an- 
ode. With the definition of the z-direction as shown in 
Fig. [1] the emission criterion becomes F tot z > 0. During 
the simulation an electron attempts to enter the system 
at a random point on the cathode where the condition 
Ftot,z > is tested by placing the electron 1 nm behind 
the surface. If true, the particle is promoted to the sur- 
face and then the test is repeated. If the result is again 
true the electron is finally emitted, with zero initial ve- 
locity. Unsuccessful emission attempts are registered and 
a series of 100 unsuccessful attempts will stop the emis- 
sion process. The inspection of F to t,z behind the surface 
ensures that all electrons already placed on the cathode 
are observed. We note that by varying the depth of the 
initial test the surface-density is affected. This depth 
can be viewed as being related to the skin-depth of the 
cathode material. 



In earlier experimental investigations by Valfells et 
al. [9] space-charge effects have been examined for photo- 
electrons extracted from a Pierce gun using a laser pulse 
with a Gaussian-like intensity profile in time. In our sim- 
ulations a radial electric field, E r = 2Vr/d 2 , represents 
the confinement of the Pierce gun, where r is the radial 
distance to the center of cathode and d equals the gap 
length that is kept constant at 2 /im. The intensity of 
a monochromatic photon pulse can be associated to the 
time interval between consecutive emission attempts, and 
a Gaussian time-interval profile can thus be applied to 
reproduce the Gaussian intensity profile. The pulse du- 
ration is defined as the time interval between first and 
last emission attempt and it is kept constant at 1 ps, 
whereas the total number of photons in a pulse (inten- 
sity) is a variable parameter. For the short-pulse limit 
under investigation, the transit time of the electrons be- 
tween the cathode and the anode exceeds the pulse du- 
ration, tp U i se << Uransit- Two series of pulse simulations 
are carried out where either the intensity of the photon 
pulse or the applied external potential is varied. When 
the intensity is varied, three regimes are simulated while 
the driving and confinement potential are kept constant 
at 2 V. For these parameters it is found that a virtual 
cathode is formed if the photon pulse exceeds a total of 
100 photons. The nomenclature we use, is to call a pulse 
with 100 photons a low intensity pulse. All higher inten- 
sities will result in yields of electrons less than 100%. For 
the medium intensity the number of photons in the pulse 
is increased by one order of magnitude and for the high 
intensity by two orders of magnitude. The main effect of 
the space-charge for a photon pulse of low intensity is a 
stretching of the electron pulse that is expected due to 
the Borsch effect. It should be noted that only for the 
low intensity do the peak value of the photon pulse and 
electron pulse coincide. For the medium and high photon 
pulse it is only the initial part that causes emission. The 
reason for this is that space-charge accumulation in the 
gap blocks subsequent electrons from being emitted, and 
the obtained electron pulse becomes a truncated Gaus- 
sian profile. Significant coupling in the z-direction and 
momentum transfer causes the electron pulse to spread 
out in both time and in velocity space when the medium 
and high intensity photon pulse are applied. A charac- 
teristic of the mentioned effects is that the electron pulse 
develops into a doubled peaked structure. 

For the simulations where the applied voltage is var- 
ied, a fixed pulse size of 350 photons is used. The applied 
potential is reduced from 10.0 V, labeled high in Fig. [2] 
by an order of magnitude, labeled medium, and finally 
by two orders of magnitude, labeled low. The confine- 
ment potential is scaled accordingly. In the regime la- 
beled high the external field dominates, whereas it be- 
comes comparable to the internal field in the medium 
case, and for the low case the internal field is dominant. 
The high potential conserves the structure of the exciting 
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FIG. 3: s Electron intensity as a function of time when a 
voltage of 1 V is applied over a 2 /im gap. A continuous pho- 
ton beam transforms into a pulsed signal of electrons with a 
frequency of 720 GHz. An absorbed electron is attributed a 
Gaussian profile in time and for the black curve (grey curve) 
a — is 800 fs (400 fs). Note that the frequency remains unaf- 
fected of profile-width and that the individual peaks are well 
separated. 

photon pulse when the electron pulse reaches the anode, 
although the electron pulse becomes slightly stretched. 
Also, only by applying the high potential is a 100% yield 
obtained. As can be seen in Fig.[2]do space-charge effects 
become significant when the applied potential decreases 
to the medium regime. At this value the electron pulse 
transforms into a double peaked structure and the yield 
decreases to 57%. As the applied potential decreases fur- 
ther the Coulomb interactions becomes dominant. A re- 
sult of this is that the initial shape of the photon pulse 
gets completely washed out and it widens fourfold. The 
yield is reduced to 24%. 

In conclusion, for both series of simulations the 
medium profiles are double-peaked structures. This char- 
acteristic of the electron pulse is consistent with experi- 
mental observations by Valfells et al. [9] where a distinct 
feature of the electron pulse in the regime where space- 
charge effects are significant is indeed a double-peaked 
structure. 

Next we look at the simulations where the cathode is 
subject to a continuous photon beam, t transit << tb eam . 
In these simulations the intensity of the photon beam 
is sufficient for the cathode to reach saturation. The 
condition for saturation is that the field created by the 
space-charge E sc , z exceeds the applied field everywhere 
on the cathode, E z < —E scz . As long as this situation 
holds the emission process is inhibited. 

Simulations of a continuous photon beam exciting the 
cathode are conducted without any confining potential in 
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FIG. 4: Frequency as a function of applied electric field. 
The deviating behavior at small gap-lengths (large electric 
fields) arises when the transit time of a charge-sheet becomes 
comparable to the period of the density oscillation. Solid line 
is a linear fit to the obtained data. Dashed lines are from the 
proposed model given by Equ. [2] 



the xy-plane. The parameters varied are either the gap 
length d or the applied voltage V. In both cases a surpris- 
ing effect is that the continuous photon beam generates a 
train of electron pulses with a well defined seperation in 
time corresponding to an oscillating current. These os- 
cillations will be the main focus in the remainder of this 
paper. Representative results for the pulsed current at 
the cathode, jc{t), and at the anode, ja(*), are shown in 
Fig. [3] All simulations with the applied electric field ex- 
ceeding 0.1 Y / fxm show clear oscillations in the electron 
density as a function of time at the cathode. At the an- 
ode the internal repulsion between electrons has caused 
an averaging wash-out dependent on the transit time. It 
should be noted that the period of the oscillations is un- 
affected by this wash-out, but a significant decrease in 
the amplitude occurs. 

In the first series of simulations the gap length is varied 
from 0.5 /^m to 4.0 /im and a constant potential of 1 V is 
applied. Only the first density peak expands freely and 
the number of electrons in this initial peak is larger than 
for the following peaks. For the following peaks the num- 
ber of electrons per pulse only shows small fluctuations. 
This latter observation indicates that the system quickly 
reaches a steady state. In Fig. [4] one should note that 
the determined frequency bends off when the gap width 
becomes small, increasing the applied field E. The other 
series of simulations conducted varies the applied voltage 
from 0.2 V to 4.0 V while keeping the gap width fixed 
at 2 /im. The correlation between the frequency of the 
density peaks and the applied field is shown in Fig. [4] 

To understand the density oscillations observed in 
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FIG. 5: Map of cathode showing regions open for emission, 
black, whereas grey and white regions are blocked. The circle 
marks the rim of the cathode. Left (1A), a blockade is formed 
as electrons are added. Middle (IB), the complete blockade 
has formed. Right (2) , a reopening of the emitter occurs when 
the electron-sheet exceeds a critical distance. 



Fig. [3] three phases of the cathode are of special impor- 
tance: (1A) Creation of a charge-sheet. (IB) Blockade, 
when a full charge-sheet has formed. (2) Reopening. For 
each of these phases we determined the total or effective 
z-component of the electric field at the cathode surface, 
which we define as E to t.z = E z + E SCtZ . The results are 
shown in Fig. [5] which is a map of the regions open for 
emission. From Fig. [5]TA and Fig. [5]TB it appears that 
the area blocked by a single electron limits to its vicinity 
within a discrete region on the surface of the cathode. 
Also to be noted is that the area blocked by an electron 
is unaffected by the presence of other electrons. 

These two observations indicate that the total number 
of electrons in a charge sheet will be determined by the 
area fraction blocked by a single electron on the cathode 
surface. When the surface is fully covered, Fig. [HJ-1B, 
the emission of electrons stops and the cathode enters 
a blocked state. Regions distant from an electron are 
subject to an enhanced blockade when the newly emit- 
ted charge-sheet starts to move toward the anode. The 
enhancement arises because the blockade from an elec- 
tron in a fixed point on the cathode depends both on 
the distance and on the angle to the surface. As the dis- 
tance increases further, the field from the charge sheet 
transforms from discrete to nearly continuous. Finally, 
when the charge sheet has covered a critical distance a 
reopening of the emitter occurs, Fig. [5]-2. At this dis- 
tance E to t,z changes sign and the formation of a new 
charge-sheet starts. Because the field distribution is now 
continuous the whole area reopens almost immediately. 
The sequence of these three states reveals that the pe- 
riod of the density oscillations is determined by the time 
it takes a charge-sheet to cover a critical distance. 

There is no simple correlation between the number of 
electrons in a charge sheet, or neither the applied po- 
tential nor the gap length. The number is between [66, 
123] when the potential is varied and [66, 218] when the 
gab length is varied, however, a clear trend in both cases 
is that the number of electrons increases as the E z in- 
creases. To understand this qualitatively it is reasonable 



to assume that the building time of a charge-sheet does 
not change significantly when E z is varied. Also, the 
speed of the sheet depends mostly on the applied field. 
The extra electrons entering a charge-sheet as the applied 
potential increases is therefore a consequence of a larger 
distance traveled by the former blocking charge-sheet be- 
fore the onset of the enhanced blockade. 

With these assumptions we can derive a simplified an- 
alytical model for the density oscillations. Suppose the 
cathode releases a thin sheet of charge which is then 
accelerated by the external field E z . Assume that the 
charge sheet is a disk of radius R, like the cathode it- 
self, situated at position z in front of the cathode. The 
field created at the center of the cathode is E sc ^ z (z) = 

Eq (l - ^ Z 2 +R2 ) , where E Q = a/2e Q and a = -Ne/nR 2 
is the charge density of a sheet containing TV electrons. 
The distance at which E tot z = in the center of the disk 



is then z = R(l - E z /Eq)/ Jl-(1- E Z /E Q ) . Let us 
admit that the acceleration of the charge sheet due to the 
external field is constant, a = eE z /m, m being the elec- 
tron mass. The time until the total field vanishes, and 
the cathode is ready to build the next sheet of charge, is 
to = y /l 2mzo/eE z and thus the frequency of the charge 
pulses is / = 1/to: 
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eE z 
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E z < E Q , then / 



If the external field is weak, 

> 1/4 3/4 

~ E z . Eq. (2) predicts that for a large 

E z -> E , the 



requency diverges, or the 



eE z 
2mR 

external field, 

oscillation period tends to zero. In that case the Coulomb 
field is significantly weaker relatively to the external field 
and the cathode is never blocked. 
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